Development and evaluation of an ozone deposition scheme for coupling to a terrestrial biosphere model

Ozone (O3) is a toxic air pollutant that can damage plant leaves and substantially affect the plant's gross primary production (GPP) and health. Realistic estimates of the effects of tropospheric anthropogenic O3 on GPP are thus potentially important to assess the strength of the terrestrial biosphere as a carbon sink. To better understand the impact of ozone damage on the terrestrial carbon cycle, we developed a module to estimate O3 uptake and damage of plants for a state-of-the-art global terrestrial biosphere model called OCN. Our approach accounts for ozone damage by calculating (a) O3 transport from 45 m height to leaf level, (b) O3 flux into the leaf, and (c) ozone damage of photosynthesis as a function of the accumulated O3 uptake over the lifetime of a leaf. A comparison of modelled canopy conductance, GPP, and latent heat to FLUXNET data across European forest and grassland sites shows a general good performance of OCN including ozone damage. This comparison provides a good baseline on top of which ozone damage can be evaluated. In comparison to literature values, we demonstrate that the new model version produces realistic O3 surface resistances, O3 deposition velocities, and stomatal to total O3 flux ratios. A sensitivity study reveals that key metrics of the air-to-leaf O3 transport and O3 deposition, in particular the stomatal O3 uptake, are reasonably robust against uncertainty in the underlying parameterisation of the deposition scheme. Nevertheless, correctly estimating canopy conductance plays a pivotal role in the estimate of cumulative O3 uptake. We further find that accounting for stomatal and non-stomatal uptake processes substantially affects simulated plant O3 uptake and accumulation, because aerodynamic resistance and non-stomatal O3 destruction reduce the predicted leaf-level O3 concentrations. Ozone impacts on GPP and transpiration in a Europe-wide simulation indicate that tropospheric O3 impacts the regional carbon and water cycling less than expected from previous studies. This study presents a first step towards the integration of atmospheric chemistry and ecosystem dynamics modelling, which would allow for assessing the wider feedbacks between vegetation ozone uptake and tropospheric ozone burden

Does the growth response of woody plants to elevated CO2 increase with temperature? A model-oriented meta-analysis

The temperature dependence of the reaction kinetics of the Rubisco enzyme implies that, at the level of a chloroplast, the response of photosynthesis to rising atmospheric CO2 concentration (Ca) will increase with increasing air temperature. Vegetation models incorporating this interaction predict that the response of net primary productivity (NPP) to elevated CO2 (eCa) will increase with rising temperature and will be substantially larger in warm tropical forests than in cold boreal forests. We tested these model predictions against evidence from eCa experiments by carrying out two meta-analyses. Firstly, we tested for an interaction effect on growth responses in factorial eCa × temperature experiments. This analysis showed a positive, but nonsignificant interaction effect (95% CI for above-ground biomass response = −0.8, 18.0%) between eCa and temperature. Secondly, we tested field-based eCa experiments on woody plants across the globe for a relationship between the eCa effect on plant biomass and mean annual temperature (MAT). This second analysis showed a positive but nonsignificant correlation between the eCa response and MAT. The magnitude of the interactions between CO2 and temperature found in both meta-analyses were consistent with model predictions, even though both analyses gave nonsignificant results. Thus, we conclude that it is not possible to distinguish between the competing hypotheses of no interaction vs. an interaction based on Rubisco kinetics from the available experimental database. Experiments in a wider range of temperature zones are required. Until such experimental data are available, model predictions should aim to incorporate uncertainty about this interaction

Vertically divergent responses of SOC decomposition to soil moisture in a changing climate

The role of soil moisture for organic matter decomposition rates remains poorly understood and underrepresented in Earth System Models (ESMs). We apply the Dual Arrhenius Michaelis-Menten (DAMM) model to a selection of ESM soil temperature and moisture outputs to investigate their effects on decomposition rates, at different soil depths, for a historical period and a future climate period. Our key finding is that the inclusion of soil moisture controls has diverging effects on both the speed and direction of projected decomposition rates (up to ± 20%), compared to a temperature-only approach. In the top soil, the majority of these changes is driven by substrate availability. In deeper soil layers, oxygen availability plays a relatively stronger role. Owing to these different moisture controls along the soil depth, our study highlights the need for depth-resolved inclusion of soil moisture effects on decomposition rates within ESMs. This is particularly important for C-rich soils in regions which may be subject to strong future warming and vertically opposing moisture changes, such as the peat soils at northern high latitudes.Vertically divergent responses of SOC decomposition to soil moisture in a changing climatepublishedVersio

Modelling sun-induced fluorescence and photosynthesis with a land surface model at local and regional scales in northern Europe

Recent satellite observations of sun-induced chlorophyll fluorescence (SIF) are thought to provide a large-scale proxy for gross primary production (GPP), thus providing a new way to assess the performance of land surface models (LSMs). In this study, we assessed how well SIF is able to predict GPP in the Fenno-Scandinavian region and what potential limitations for its application exist. We implemented a SIF model into the JSBACH LSM and used active leaf-level chlorophyll fluorescence measurements (Chl F) to evaluate the performance of the SIF module at a coniferous forest at Hyytiala, Finland. We also compared simulated GPP and SIF at four Finnish micrometeorological flux measurement sites to observed GPP as well as to satellite-observed SIF. Finally, we conducted a regional model simulation for the Fenno-Scandinavian region with JSBACH and compared the results to SIF retrievals from the GOME-2 (Global Ozone Monitoring Experiment-2) space-borne spectrometer and to observation-based regional GPP estimates. Both observations and simulations revealed that SIF can be used to estimate GPP at both site and regional scales. At regional scale the model was able to simulate observed SIF averaged over 5 years with r(2) of 0.86. The GOME-2-based SIF was a better proxy for GPP than the remotely sensed fA-PAR (fraction of absorbed photosynthetic active radiation by vegetation). The observed SIF captured the seasonality of the photosynthesis at site scale and showed feasibility for use in improving of model seasonality at site and regional scale.Peer reviewe

A few extreme events dominate global interannual variability in gross primary production

Understanding the impacts of climate extremes on the carbon cycle is important for quantifying the carbon-cycle climate feedback and highly relevant to climate change assessments. Climate extremes and fires can have severe regional effects, but a spatially explicit global impact assessment is still lacking. Here, we directly quantify spatiotemporal contiguous extreme anomalies in four global data sets of gross primary production (GPP) over the last 30 years. We find that positive and negative GPP extremes occurring on 7% of the spatiotemporal domain explain 78% of the global interannual variation in GPP and a significant fraction of variation in the net carbon flux. The largest thousand negative GPP extremes during 1982–2011 (4.3% of the data) account for a decrease in photosynthetic carbon uptake of about 3.5 Pg C yr−1, with most events being attributable to water scarcity. The results imply that it is essential to understand the nature and causes of extremes to understand current and future GPP variability

Predicting resilience through the lens of competing adjustments to vegetation function

There is a pressing need to better understand ecosystem resilience to droughts and heatwaves. Eco-evolutionary optimization approaches have been proposed as means to build this understanding in land surface models and improve their predictive capability, but competing approaches are yet to be tested together. Here, we coupled approaches that optimize canopy gas exchange and leaf nitrogen investment, respectively, extending both approaches to account for hydraulic impairment. We assessed model predictions using observations from a native Eucalyptus woodland that experienced repeated droughts and heatwaves between 2013 and 2020, whilst exposed to an elevated [CO2] treatment. Our combined approaches improved predictions of transpiration and enhanced the simulated magnitude of the CO2 fertilization effect on gross primary productivity. The competing approaches also worked consistently along axes of change in soil moisture, leaf area, and [CO2]. Despite predictions of a significant percentage loss of hydraulic conductivity due to embolism (PLC) in 2013, 2014, 2016, and 2017 (99th percentile PLC > 45%), simulated hydraulic legacy effects were small and short-lived (2 months). Our analysis suggests that leaf shedding and/or suppressed foliage growth formed a strategy to mitigate drought risk. Accounting for foliage responses to water availability has the potential to improve model predictions of ecosystem resilience

On the use of Earth Observation to support estimates of national greenhouse gas emissions and sinks for the Global stocktake process: lessons learned from ESA-CCI RECCAP2

The Global Stocktake (GST), implemented by the Paris Agreement, requires rapid developments in the capabilities to quantify annual greenhouse gas (GHG) emissions and removals consistently from the global to the national scale and improvements to national GHG inventories. In particular, new capabilities are needed for accurate attribution of sources and sinks and their trends to natural and anthropogenic processes. On the one hand, this is still a major challenge as national GHG inventories follow globally harmonized methodologies based on the guidelines established by the Intergovernmental Panel on Climate Change, but these can be implemented differently for individual countries. Moreover, in many countries the capability to systematically produce detailed and annually updated GHG inventories is still lacking. On the other hand, spatially-explicit datasets quantifying sources and sinks of carbon dioxide, methane and nitrous oxide emissions from Earth Observations (EO) are still limited by many sources of uncertainty. While national GHG inventories follow diverse methodologies depending on the availability of activity data in the different countries, the proposed comparison with EO-based estimates can help improve our understanding of the comparability of the estimates published by the different countries. Indeed, EO networks and satellite platforms have seen a massive expansion in the past decade, now covering a wide range of essential climate variables and offering high potential to improve the quantification of global and regional GHG budgets and advance process understanding. Yet, there is no EO data that quantifies greenhouse gas fluxes directly, rather there are observations of variables or proxies that can be transformed into fluxes using models. Here, we report results and lessons from the ESA-CCI RECCAP2 project, whose goal was to engage with National Inventory Agencies to improve understanding about the methods used by each community to estimate sources and sinks of GHGs and to evaluate the potential for satellite and in-situ EO to improve national GHG estimates. Based on this dialogue and recent studies, we discuss the potential of EO approaches to provide estimates of GHG budgets that can be compared with those of national GHG inventories. We outline a roadmap for implementation of an EO carbon-monitoring program that can contribute to the Paris Agreement